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INTRODUCTION: 


This  study  involves  the  technical  development  and  clinical  testing  of  a  novel  technique  for 
magnetic  resonance  imaging  (MRI)  guided  prostate  biopsy  in  a  1.5T  horizontal  bore  scanner 
using  a  dedicated  interventional  table.  We  primarily  hypothesized  that  the  integration  of 
diagnostic  and  interventional  MRI  enables  needle  biopsy  targeting  to  foci  of  tumor  recurrence 
after  radiotherapy,  and  will  enable  a  determination  of  the  diagnostic  accuracy  of  MRI  in  mapping 
sub-sites  of  tumor  recurrence  after  radiotherapy.  In  this  final  report,  we  provide  strong  evidence 
supporting  our  primary  study  hypothesis. 

The  clinical  trial  was  designed  to  enroll  patients  with  suspicion  of  locally  recurrent  prostate  cancer  after 
external  beam  radiotherapy  In  a  pilot,  three-stage  trial  design.  The  first  stage  included  technical 
development  and  optimization  of  the  Imaging  technique.  The  second  stage  included  technical 
development  and  optimization  of  the  biopsy  technique,  followed  by  a  clinical  evaluation  stage.  Stages  1 
and  2  have  completed  accrual,  where  16  patients  were  enrolled  In  stage  1 ,  and  4  patients  were  enrolled 
in  stage  2.  Stage  3  continues  to  accrue,  and  to  date  has  enrolled  14  patients  of  a  planned  target  accrual 
of  30  patients.  A  companion  pilot  study  investigating  profiles  of  hypoxia  has  completed  accrual  with  13 
patients  and  is  now  closed.  The  clinical  evaluation  stage  3  of  the  trial  is  expected  to  close  to  accrual 
within  the  next  12  months. 

BODY: 

USAMRCMC  approval  for  the  conduct  of  the  clinical  trial  was  obtained  on  October  13,  2006,  therefore  the 
statement  of  work  was  delayed  by  one  year  for  regulatory  approval.  Here,  we  report  final  and  most 
updated  results  for  each  task  of  the  project. 

Aim  1:  To  test  the  hypothesis  that  stereotactic  MRI-guided  prostate  needle  biopsies  can  be  performed 
with  improved  anatomic  targeting  accuracy  and  tolerability  when  patients  are  positioned  supine  in  a 
conventional  MR  scanner. 

Task  la.  To  design  a  custom  interventional  MRI  table  and  hardware  system  that  provides  immobilization 
and  perineal  access  in  the  supine  position. 

The  table  and  hardware  system  has  been  designed  in  collaboration  with  industry  (Sentinelle  Medical 
Inc.),  and  our  current  system  (depicted  below)  is  the  third  prototype  now  scheduled  for  product  release  in 
2010.  The  navigation  system  provides  access  to  the  perineum  in  an  immobilized  supine  patient  position 
within  a  1 .5T  GE  Signa  scanner,  and  a  stereotactic  transperineal  template  assembly  which  Integrates  a 
commercial  endorectal  Imaging  coil  (C  -  Medrad  Inc.).  The  pelvis  and  legs  are  immobilized  using  velcro 
straps. 


Task  1b.  To  measure  prostate  motion  while  patients  are  immobiiized  on  the  custom  table. 

Prostate  motion  was  measured  by  identifying  3D  fiducial  reference  points  visible  on  both  the  registration 
T2  images  at  the  beginning  of  the  procedure,  and  again  on  the  first  and  last  needle  verification  image. 

The  majority  of  reference  points  consisted  of  implanted  gold  marker  visualized  as  an  area  of  hypointensity 
on  T2  FSE  images.  Where  gold  markers  were  not  present,  distinct  anatomic  points  of  hyper  or  hypo- 
intensity  were  used.  Prostate  motion  was  determined  by  averaging  displacements  of  all  visible  points  in  a 
given  image  pair.  Between  the  T2  images  used  for  guidance  and  needle  verification  images  after 
insertion  of  the  first  biopsy  needles,  median  vector  displacement  for  the  entire  cohort  was  1 .3  mm  (SD 

I . 0,  Max  3.2).  The  largest  displacements  were  observed  in  the  SI  dimension,  where  image  resolution 
was  3mm,  and  along  the  trajectory  of  force  applied  by  needles.  Similarly,  there  was  minimal  prostate 
displacement  between  the  registration  T2  images  and  final  needle  verification  image,  with  a  median  (SD) 
vector  displacement  of  3.1  mm  (1 .5),  within  Image  spatial  resolution.  This  data  support  the  hypothesis 
that  our  positioning  and  sedation  protocol  minimizes  prostate  motion,  and  enables  stereotactic  needle 
targeting  under  the  assumption  of  a  static  system. 

Task  1c.  To  measure  tolerance  and  feasibility  of  the  integrated  diagnostic  and  interventional  procedure. 

Patient  in  the  evaluation  cohort  (stage  3)  underwent  the  integrated  diagnostic  and  biopsy  procedure 
under  propofol  sedation.  As  such,  pain  and  tolerability  VAS  scores  were  reported  to  be  0  in  all  but  one 
patient,  who  reported  a  0.5  score  (out  of  10  points,  0=no  pain/discomfort,  10=worst  possible 
pain/discomfort).  Regarding  adverse  events,  one  patient  experienced  a  Grade  3  pulmonary  embolism 
within  a  few  days  of  the  biopsy  procedure.  He  was  investigated  and  found  to  have  a  hypercoagulative 
syndrome,  and  the  event  was  deemed  unlikely  to  be  related  to  the  procedure.  A  second  patient 
experienced  a  Grade  2  vasovagal  hypotensive  event  1  hour  after  the  procedure  which  resolved  with 
conservative  measures.  No  patients  experienced  an  infection  or  other  complications  related  to  the  biopsy 
procedure. 

Feasibility  metrics  included  the  following:  frequency  of  visible  suspicious  tumor  targets  for  biopsy, 
frequency  of  complete  biopsy  sampling,  and  procedure  times.  Of  14  patients  enrolled  to  date  in  stage  3 
of  the  trial,  one  procedure  was  aborted  due  to  MRI  scanner  malfunction.  All  thirteen  evaluable  patients 
had  regional  MRI  features  suspicious  for  cancer  and  specifically  targeted  for  biopsy.  Intended  biopsy 
sampling  (tumor  targets  +  a  minimum  of  sextant  sampling)  was  successfully  achieved  in  10  patients. 
Insufficient  diagnostic  sampling  was  observed  in  the  first  two  patients  due  to  poor  core  quality,  which  led 
to  a  change  manual  handling  of  the  MRI  compatible  biopsy  gun.  A  later  patient  had  insufficient  sampling 
due  longer  needle  length  requirements  (thick  perineal  fat),  and  subsequent  procedures  utilized  longer 
biopsy  guns  (175mm  instead  of  150mm).  The  number  of  biopsy  cores  sampled  per  patient  ranged  from 
8-13.  Mean  diagnostic  MRI  +  biopsy  procedure  time  was  76  min  (56-93  min).  For  each  needle  verification 
Interval,  two  biopsies  were  acquired  from  targets  in  the  right  and  left  lobes  of  the  gland.  The  mean  time 
for  needle  Insertion,  needle  verification,  and  biopsy  acquisition  fora  biopsy  interval  was  7.5  min  (5.8- 

I I . 6min).  After  patient  5,  set-up  was  performed  in  the  MRI  scanner  suite.  Trends  in  time  performance 
are  depicted  below  through  the  course  of  the  study. 


Task  1d.  To  confirm  the  geometric  needle  targeting  accuracy  of  the  supine  trans-perineal  system  in  a  gei 
phantom. 

A  phantom  experiment  was  conducted  using  Aegis  targeting  software,  and  Prostate  Phantom  (Modei  053, 
CiRS  inc.),  as  well  as  Invivo  16  gage  automatic  biopsy  gun.  Mean  needle  targeting  accuracy  (n=7)  was: 
x=0.34  +/-  0.31  mm,  and  y=-1 .57+/-1 .09  mm.  A  posterior  systematic  but  submillimeter  bias  can  be 
addressed  by  re-caiibrating  the  device.  This  phantom  data  was  reported  at  the  ISMRM  Annual  meeting  in 
2007.  (1,2)  Resuits  indicate  high  performance  of  the  needie  navigation  system,  and  was  deemed 
sufficiently  accurate  to  proceed  with  clinical  testing. 

Task  1e.  To  test  the  geometric  and  anatomic  needie  targeting  accuracy  of  the  supine  trans-perineai 
system  in  patients  with  suspected  iocal  persistence  of  prostate  cancer  after  external  beam  radiotherapy. 

The  navigation  system  is  a  fundamental  component  of  image  guided  prostate  cancer  targeting.  The 
following  figure  illustrates  the  feedback  relationships  designed  to  improve  prostate  cancer  targeting 
accuracy  of  the  navigation  system. 
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Figure  1  ■  Illustration  of  the  feedback  relationships  designed  to  improve  targeting  accuracy  of  the  navigation 
system.  General  MRI  features  of  prostate  cancer  (nodular  area  of  hypointensity  on  T2  weighted  images,  low 
apparent  diffusion  coefficient,  elevated  choline  relative  to  citrate,  and  intense  and  rapid  contrast 
enhancement  and  washout  visualized  here  in  the  right  peripheral  zone)  are  used  to  identify  suspicious 
Imaging  targets  for  biopsy  [red  target).  30  needle  verification  images  are  used  to  assess  and  improve  the 
stereotactic  performance  of  the  navigation  system  by  comparing  the  actual  needle  location  (dark  signal  void) 
to  that  intended  (red  target]  in  stereotactic  space.  These  images,  in  conjunction  with  deformable  image 
registration,  can  also  provide  metrics  of  anatomic  needle  targeting  accuracy.  This  data  can  be  used  to  create 
needle  deformation  models  that  could  be  tncorporated  into  the  navigation  system  for  rmproved  needie 
targeting.  Most  importantly,  3D  needle/tissue  verifFcation  images  and  deformable  image  registration  enable 
accurate  voxel-to-tissueco-iocalizadon  that  will  validate  and  help  refine  quantitative  MRI  parameters 
diagnostic  of  cancer  burden  and  bioiogical  virulence,  thereby  improving  imaging  target  accuracy 


The  inplane  stereotactic  targeting  accuracy  was  calculated  by  finding  the  distance  between  the  actual 
needle  and  the  intended  target  location  in  MRI  coordinates.  The  actual  needle  location  was  recorded  as 
the  center  of  the  signal  void  visible  in  the  needle  verification  images.  For  a  total  of  121  needle  targets  in 
12  patients,  the  mean  (SD)  in  plane  targeting  accuracy  in  millimeter  was  x=0.44  (1 .77),  y=1 .46  (1 .77),  and 
vector  2.57  (1 .41 ).  The  inplane  error  was  distributed  randomly  in  the  right  to  left  direction  but  a  small 
posterior  systematic  bias  of  the  inplane  errors  of  approximately  1  mm  was  observed.  The  following  figure 
depicts  all  stereotactic  needle  target  errors  as  well  as  the  mean  (SD)  for  the  entire  cohort,  with  each 
needle  considered  an  independent  variable.  As  needle  targets  within  a  given  patient  are  unlikely  to  be 
independent  variables,  we  also  computed  the  mean  (SD)  for  each  patient.  Across  patients,  the  mean 
vector  targeting  error  was  2.5  mm  (range  1 .2  -  4.7  mm). 
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In  keeping  with  prostate  motion  data,  hardware  motion  was  minimal  throughout  the  biopsy  procedures, 
with  a  mean  hardware  motion  of  0.64mm  (SD  0.39mm),  lending  credence  to  the  stereotactic  guidance 
assumptions  of  a  stable  system. 


Stereotactic  Device  Motion 
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Anatomic  targeting  accuracy  will  be  determined  by  measuring  the  vector  discrepancy  between  the 
intended  anatomic  target  for  biopsy,  and  sampled  region  as  visualized  on  needle  verification  images. 
Analysis  is  ongoing,  using  3D  maps  and  FEM  deformable  registration  techniques  based  on  the  prostate 
surface  boundary  manually  delineated  on  MRI. 

It  is  our  impression  based  on  this  data  that  the  majority  of  targeting  inaccuracies  occur  due  to  needle 
deflection  at  the  prostate  boundary,  which  is  more  pronounced  for  targets  located  at  the  margins  of  the 
prostate  gland.  Ongoing  analysis  will  attempt  to  model  and  predict  for  needle  deflection  in  order  to 
improve  overall  targeting  accuracy  of  the  system. 

Gland  deformation  is  minimal  with  a  transperineal  approach,  but  needle  insertion  does  slightly  translate 
the  gland  superiorly,  as  previously  reported  by  others.  Further  work  will  continue  to  refine  accurate  and 
responsive  registration  of  the  stereotactic  system,  model  and  account  for  needle  deflection  and  prostate 
translation  for  improved  targeting  performance. 


Aim  2:  To  test  the  hypothesis  that  stereotactic  MRI-guided  needle  biopsies  will  enable  the  spatial 
delineation  of  local  prostate  cancer  persistence  after  radiotherapy. 

Task  2a.  To  create  3D  maps  of  malignant  biopsy  sites  in  reference  to  MRI  prostate  anatomy. 

3D  maps  of  malignant  biopsy  sites  in  reference  to  MRI  prostate  anatomy  have  been  successfully 
generated  in  all  patients.  A  case  example  is  depicted  in  the  following  figure.  In  b)  navigation  software 
(Aegis,  Sentinelle  Medical  Inc.)  displays  needle  target  (green)  onto  diagnostic  T2  images.  In  c)  needle 
verification  images  document  the  actual  location  of  3  biopsy  needles  (signal  voids  R  and  L)  in  reference 
to  intended  target  (green).  In  d),  a  map  of  benign  (green),  malignant-dense  (red),  malignant-microfocus 
(pink),  and  indeterminate  (white)  biopsy  cores  are  depicted  in  reference  to  imaging  features  suspicious  for 
tumor  on  DCE  (light  blue),  ADC  (purple),  and  T2  (royal  blue).  In  this  case,  the  boundary  of  the  T2-visible 
tumor  was  later  revised  posteriorly  based  on  a  benign  biopsy  core  in  this  region. 


This  tumor  map  can  be  rendered  in  3D,  as  depicture  below  for  a  different  patient  with  a  tumor  visualized 
and  pathologically  confirmed  in  the  right  lateral  lobe  of  the  gland. 


Task  2b.  To  correlate  biopsy  tissue  histology  with  spatially  corresponding  diagnostic  MR!  data. 

The  median  age  at  study  accrual  was  72  years  (63-84).  All  received  RT  (+/-  hormones)  for  intermediate 
(n=13)  or  low-risk  (n=2)  prostate  cancer  at  initial  diagnosis.  The  median  EBRT  dose  was  75.6Gy  (60-79.8 
Gy).  Only  2  patients  received  neoadjuvant  hormone  therapy  for  a  period  of  5  months.  The  median  PSA  at 
the  time  of  biopsy  was  3.8  ng/ml  (2. 9-7. 9).  The  median  time  from  RT  to  biopsy  was  5.8  years  (1 .7-1 1 .4). 

Of  13  evaluable  patients  evaluable  to  date,  10  patients  were  found  to  have  recurrent  cancer  on  their  MRI- 
guided  biopsy.  The  majority  of  patients  had  Gleason  score  7  recurrence  disease  (n=7,  Gleason  score 
range  6-9).  The  site  of  recurrence  corresponded  to  suspicious  imaging  findings  in  all  but  2  patients.  In 
one  patient  the  malignant  biopsy  did  not  spatially  correspond  with  a  suspicious  imaging  focus.  In  the 
second  patient,  disease  was  diffusely  distributed  throughout  the  gland  and  no  focus  of  disease  could  be 
delineated.  In  all  patients  with  focal  recurrences,  only  1  focus  of  disease  was  present.  Extracapsular 
extension  of  disease  was  identified  in  2  patients. 

Ongoing  analysis  will  have  2  separate  expert  observers  delineate  regions  that  they  consider  suspicious 
for  tumor-burden  based  on  multiparametric  MRI.  Overall  accuracy  (sensitivity/specificity)  of  the  experts 
will  be  determined  based  on  co-registered  biopsy  data.  (Collaboration  M.  Haider) 

Task  2c.  To  define  a  profile  of  MRI  measurements  that  accurately  delineates  sites  of  tumor  persistence 
within  the  prostate  gland  after  external  beam  radiotherapy. 

Imaging  of  actual  needle  core  locations  enabled  the  identification  of  corresponding  imaging  voxels  on 
diagnostic  images.  Once  the  corresponding  3D  location  was  identified  on  T2,  a  1x1  mm  area 
(corresponding  to  biopsy  size  and  representing  1 .2  voxels)  was  delineated  on  sequential  slices 
corresponding  to  the  reported  length  of  biopsy  core.  The  average  and  standard  deviation  of  voxel  values 
contained  within  each  volume  was  calculated.  Mean  T2  was  then  derived  for  each  biopsy  location.  While 
up  to  3  slices  (representing  the  entire  biopsy  length)  were  included  to  derive  quantitative  metrics  for 
benign  cores,  only  slices  corresponding  to  tumor  (as  estimated  from  distance  between  inked  end  to 
tumor)  were  used  to  derive  quantitative  metrics  for  the  malignant  cores.  Quantitative  metrics  derived  from 
biopsies  containing  microfoci  (defined  as  a  single  focus  less  than  or  equal  to  a  40  x  microscopic  field) 
(Allan  RW,  Journal  of  Urology,  2003;)  or  less  than  5%  of  involved  core  were  excluded  from  analysis. 

Six  patients  with  quantitative  T2  images  acquired  at  the  time  of  biopsy  have  been  analyzed  to  date,  for  a 
total  of  45  sampled  regions.  There  was  a  statistically  significant  difference  in  mean  T2  between  benign 
(85ms)  and  malignant  (1 14ms)  cores  (p=0.02),  although  substantial  overlap  remains  in  the  range  of  T2 
values,  and  no  threshold  value  can  accurately  distinguish  benign  and  malignant  disease  (following  figure). 
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A  similar  quantitative  analysis  is  ongoing  to  investigate  the  potential  for  other  quantitative  metrics  (Ktrans, 
ADC,  Choline+Creatine/Citrate)  to  distinguish  benign  and  malignant  tissue  voxels.  Multivariate  analysis 
will  then  proceed  to  determine  which  combination  of  quatitative  threshold  achieves  the  highest  diagnostic 
accuracy  in  delineating  tumor. 

Aim  3:  To  test  the  hypothesis  that  integrated  diagnostic  and  interventional  MRI  will  enable  the  biological 
characterization  of  prostate  cancer  persistence  within  the  prostate  gland  after  radiotherapy. 

Task  3a.  To  obtain  stereotactic  measures  of  tumor  hypoxia. 

Thirteen  patients  received  an  infusion  of  Pimonidazole  (0.5mg/m^)  over  20  minutes  approximately  18 
hours  prior  to  the  biopsy  to  serve  as  a  immunohistochemical  biomarker  of  tumor  hypoxia.  An  example  of 
a  needle  core  acquired  under  stereotactic  MRI  guidance  and  stained  with  an  antibody  to  Pimonidazole  is 
depicted  below.  Staining  and  analysis  of  the  complete  cohort  in  ongoing  (collaboration  Dr.  R.  Bristow). 


T2*/Blood  Oxygen  Level  Dependent  (BOLD)  MRI  image  sequences  are  thought  to  provide  estimates  of 
tissue  oxygenation  depending  on  variable  magnetic  properties  of  oxy  and  deoxyhemoglobin.  All  13 
patients  in  the  companion  cohort  had  quantitative  T2*  images  acquired  immediately  prior  to  the  biopsy 
(image  example  below).  One  Pimonidazole  staining  and  quantification  is  complete,  voxel/pathological 
correlation  will  proceed  using  the  methodology  described  in  Task  2c. 


KEY  RESEARCH  ACCOMPLISHMENTS: 

•  Design  of  a  system  for  integrated  diagnostic  and  interventional  MRI  using  stereotactic 
transperineal  needle  navigation. 

•  Clinical  demonstration  of  feasibility,  tolerance,  target  stability,  and  targeting  accuracy  of  the 
procedure. 

•  Raising  the  importance  and  value  of  3D  verification  imaging  in  imaging/pathological  correlative 
investigations. 

•  Demonstration  that  recurrent  prostate  cancer  after  radiotherapy  can  usually  be  localized  to  a 
single  focus  of  disease  visible  on  MRI,  which  can  be  targeted  for  biopsy. 

•  Preliminary  evidence  that  quantitative  T2  may  contribute  in  distinguishing  tumor  from  benign 
tissue. 
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CONCLUSIONS: 

Feasibility  of  integrating  multi-parametric  MRI  for  online  guidance  of  prostate  biopsy  is  demonstrated. 

The  value  of  3D  imaging  to  document  actual  location  of  biopsy  cores  in  reference  to  anatomic  boundaries 
cannot  be  overstated.  Online  MRI  needle  guidance  systems  with  accurate  and  responsive  navigation  help 
better  define  cancer  features  on  MRI  and  enable  tumor-targeted  diagnostics  and  therapeutics.  Future 
work  will  complete  the  ongoing  analysis,  and  expand  investigations  to  a  broader  cohort  of  patients. 
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A  System  for  Prostate  Intervention  in  a  1.5  T  MRI  Scanner  in  the  Supine  Position 


Bootsma  GJ,  Krieger  A,  lordachita  II,  Piron  C,  Richmond  J,  Sela  G,  Filleti  M,  Sie  F,  Rocca 
C,  Kirilova  A,  Brock  K,  Jaffray  DA,  Haider  MA,  Menard  C 

Introduction:  Prostate  cancer  is  the  most  likely  cancer  to  develop  in  a  Canadian  male  with  a  prevalence 
of  1  in  7.1  cases  in  2005  [1],  In  the  United  States  230,000  new  cases  of  prostrate  cancer  have  been 
projected  for  2005  [2],  Despite  improvements  in  the  delivery  and  reduction  in  associated  toxicity  of 
external  beam  radiotherapy,  local  persistence  or  recurrence  of  disease  remains  prevalent  in  25-51  %  of 
patients  [3,4],  A  careful  and  thorough  investigation  of  the  spatial  distribution  of  cancer  within  the  prostate 
gland  is  paramount  to  meaningful  progress  in  effective  prospective  management  and  therapy  for  salvage. 
Magnetic  Resonance  Imaging  (MRI)  provides  the  ability  to  deliver  excellent  soft  tissue  contrast  and 
resolution  of  the  anatomy.  It  also  provides  an  opportunity  to  spatially  characterize  pathology  and  biology 
through  dynamic  contrast-enhanced  (DCE)  imaging,  MR  spectroscopic  imaging  of  tissue  metabolites,  and 
diffusion-weighted  (DWI)  imaging.  Previous  systems  have  been  developed  for  prostrate  interventions  in  a 
standard  1.5  T  MRI  scanner  with  excellent  biopsy-needle  targeting  accuracy,  but  have  required  the 
patient  to  be  placed  in  the  prone  or  left  lateral  decubitus  position  [3],  thereby  compromising  stability, 
patient  comfort,  and  safety.  A  number  of  recent  studies  have  shown  that  prostate  motion  is  greatly 
reduced  when  patients  are  placed  in  the  supine  position  due  to  greater  patient  comfort  and  reduced 
respiratory  motion  [4-6].  In  this  work,  the  feasibility  of  creating  adequate  perineal  exposure  for  prostate 
interventions  in  the  supine  position  using  dedicated  table  architecture  is  investigated. 


Figure  1:  Picture  of  the  custom  imaging  table  for  patient 
positioning  in  the  supine  position  docked  to  GE  1.5  T  MRI 
scanner. 


Figure  2:  Image  showing  the  accessibility  of  the  perineum  and 
accommodation  of  stereotactic  perineal  template  and  coil. 


Materials/Methods:  To  permit  the  patient  to  be  positioned  supine  and  allow  access  to  the  perineum  a 
dedicated  prostate  interventional  table  (Figure  1)  was  designed  in  collaboration  with  Sentinelle  Medical 
Inc.  (Toronto,  Canada).  The  table  is  a  modification  of  Sentinelle  Medical's  Vanguard  System  which  is 
used  in  breast  imaging  and  intervention.  The  table  has  been  revised  to  provide  boot  supports  for  hip  and 
knee  flexion  and  immobilization,  pelvic  immobilization,  imaging  coil  integration,  and  ample  room  for 
perineal  exposure  and  operative  space  (Figure  2).  The  table  has  also  been  designed  to  accommodate 
both  the  “access  to  prostate  tissue  under  MRI  guidance”  (APT-MRI)  device  [7,8]  as  well  as  a  stereotactic 
perineal  template  and  coil  system  [3].  The  interventional  table  provides  ease  of  mobility  between  the  MRI 
suite  and  adjacent  rooms.  Patients  can  then  be  set-up  on  the  table  prior  to  it  being  docked  to  the  MRI 
scanner,  thus  freeing  the  unit  for  standard  clinical  workflow.  An  REB/HSRBB  approved  protocol  is 
currently  enrolling  patients  to  develop  the  revised  interventional  techniques  and  measure  organ  motion 
and  patient  comfort  in  the  new  system.  Prostate  motion  was  measured  using  a  2D  FIESTA  with  fat 
saturation  (TE/TR  -  1.9/6.2ms;  matrix  320x224,  FOV  26  cm,  8mm  slice  thickness)  and  temporal 
resolution  2sec,  and  analyzed  with  manual  POI  tracking  using  MIPAV  software  (Medical  Image 
Processing  Analysis  and  Visualization  -  NIH). 


Results:  On  an  initial  trial,  the  patient  was  placed  on  the  table  and  imaged  for  a  period  of  60  minutes  with 
the  endorectal  coil  in  place.  There  was  no  report  of  discomfort  from  the  patient,  and  perineal  exposure 
was  excellent  (Figure  2).  Respiratory  and  peristaltic  motion  observed  during  the  2D  FIESTA  sequence 


caused  negligible  movement  of  the  prostate  gland.  Maximum  displacement  of  the  prostate  gland  was 
0.51mm  and  1.02mm,  associated  with  a  mean  displacement  of  0.00  +/-0.31mm  and  0.00  +/-0.31mm  in 
the  AP  and  SI  directions,  respectively. 

Conclusions:  Our  new  system  for  prostate  intervention  provides  a  means  by  which  MRI  data  can  be 
used  to  accurately  guide  intervention  as  well  as  provide  feedback  to  validate  new  imaging  techniques  with 
the  patient  placed  in  the  supine  position  minimizing  patient  discomfort  and  prostate  motion. 

References:  [1]  NCIC:  Canadian  Cancer  Statistics  2005,  Toronto,  Canada  2005.  [2]Jemal  A.  CA  Cancer 
J  Clin,  2005:55(1  ):10-30.  [3]  Susil  RC.  Magn  Reson  Med  2004;52(3);683-7.  [4]  Mah  D.  Int  J  Radiat  Oncol 
Biol  Phys  2002:54(2):568-75.  [5]  Kitamura  K.  Int  J  Radiat  Oncol  Biol  Phys  2002;53(5):  1 117-23.  [6] 

Bayley  AJ.  Radiother  Oncol  2004;70(1):37-44.  [7]  Susil  RC.  Radiology  2004;228(3):886-94.  [8]  Krieger  A. 
IEEE  Trans  Biomed  Eng  2005;  52(2):306-13. 


A  Supine  System  for  Needle-Based  Prostate  Interventions  in  a  Cylindrical  MRI  Scanner 
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Introduction:  A  careful  and  thorough  investigation  of  the  spatial  distribution  of  cancer  within  the 
prostate  gland  is  paramount  to  progress  in  effective  tumour-targeted  therapy.  MRI  provides  an 
opportunity  to  characterize  the  distribution  of  glandular  pathology  and  biology.  Previous 
interventional  MRI  systems  have  positioned  patients  in  prone  or  left  lateral  decubitus,  thereby 
compromising  stability,  patient  comfort,  and  safety.  In  the  first  aim  of  this  work,  we  sought  to 
determine  the  feasibility  of  creating  adequate  perineal  exposure  for  prostate  interventions  in  the 
supine  position  using  dedicated  table  architecture. 

Methods:  An  interventional  table  was  designed  providing  leg  and  pelvic  immobilization,  imaging 
coil  integration,  fixation  arms  for  interventional  hardware,  operative  space  with  perineal 
exposure,  and  ease  of  transfer  between  MRI  and  other  interventional  suites.  Patients  were 
enrolled  on  the  'run-in'  phase  of  a  clinical  trial  in  order  to  develop  the  interventional  techniques 
and  measure  organ  motion  in  the  new  system.  Prostate  motion  was  measured  using  sagittal  2D 
FIESTA  (temporal  resolution  2sec),  and  analyzed  using  manual  POI  tracking. 

Results:  Six  patients  with  suspicion  of  local  recurrence  after  EBRT  have  been  enrolled  and 
investigated  on  the  system  to  date.  Exposure  was  adequate  for  transperineal  intervention  in  all 
cases.  Respiratory  and  peristaltic  motion  observed  in  one  patient  during  the  2D  FIESTA 
sequence  caused  negligible  movement  of  the  prostate  gland,  with  maximum  and  mean  (SD) 
displacements  of  O.Slmm  and  0mm  (0.31mm)  in  the  AP  direction,  and  1.02mm  and  0mm 
(0.46mm)  in  the  SI  dimension.  Recurrent  tumour  was  visualized  on  MRI,  and  confirmed  on 
subsequent  TRUS-guided  biopsy  in  5/6  patients. 

Conclusions:  The  new  system  for  prostate  intervention  provides  an  improved  means  for  online 
MRI  guidance  of  needlebased  interventions  (biopsy  and  brachy therapy).  In  next  steps,  MRI- 
guided  prostate  biopsies  will  be  performed  in  patients  suspected  of  local  recurrence  after 
radiotherapy. 
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Purpose 

Despite  a  strong  rationale  and  enthusiasm  for  distinguishing  regions  of  tumor-dense  burden 
within  the  prostate  gland  for  specific  dose  delivery,  methods  and  tools  for  clinical  investigation 
have  not  yet  matured.  A  technique  is  presented  for  integrating  pathologically  referenced 
multiparametric  MRI  using  deformable  registration  for  HDR  brachytherapy  planning. 

Methods 

Four  patients  enrolled  on  a  prospective  trial  of  tumor-targeted  MRI-guided  HDR  salvage 
brachytherapy  were  evaluated,  each  receiving  two  fractions  over  10  days.  Multiparametric 
online  MRI-guided  prostate  biopsy  with  3D  image  verification  of  needle  core  locations  enabled 
cores  and  suspicious  tumor  volumes  to  be  mapped  onto  T2  MRI.  The  GTV  was  defined  as  the 
shared  boundary  of  suspicious  image  features  and  malignant  biopsy  cores.  On  treatment  day, 
radiation  oncologist  estimated  the  intended  GTV  boundary  based  on  common  anatomical  MRI 
landmarks.  Inverse  treatment  planning  utilized  IPSA  objectives  to  achieve  tumor  V 1  lGy>95% 
and  prostate  V8Gy>95%,  while  respecting  OAR  constraints. 

MORFEUS,  a  biomechanical  model-based  deformable  image  registration  technique,  was  used  to 
perform  deformable  registration  of  diagnostic  images  acquired  during  biopsy  onto  the 
brachytherapy  planning  image.  Common  3D  points  in  the  prostate  were  identified  to  determine 
the  absolute  error  between  the  observed  and  MORFEUS-predicted  point  displacements.  The 
registered  GTV  was  compared  with  the  estimated  GTV  using  Dice's  coefficient  and  tumor  D98 
and  VI 1  Gy. 

Result 

The  average  absolute  error  of  point  displacement  was  0.5mm  LR,  1.3mm  AP  and  1.3mm  SI. 

The  average  Dice's  coefficient  was  0.75.  The  majority  of  the  geometric  discrepancy  occurred  in 


the  SI  dimension  (Figure  1,  registered  and  estimated  GTV  in  blue  and  orange,  respeetively).  The 
registered  GTV  reeeived  higher  than  intended  dose  in  5/7  fraetions  (D98  range  lO.OGy  - 
13.  IGy).  Tumor  coverage  was  compromised  in  2/7  fractions  with  the  estimated  technique 
(VI  IGy  range  85.6%  -  100%). 

Conclusion 

Feasibility  of  integrating  pathologically  referenced  multi-parametric  MRI  for  GTV  delineation 
during  HDR  prostate  brachytherapy  has  been  demonstrated.  The  accuracy  achieved  with  the 
deformable  registration  was  smaller  than  the  largest  voxel  dimension.  A  small,  but  potentially 
relevant  dosimetric  discrepancies  between  estimated  and  deformably  registered  GTV  was 
observed.  Ongoing  work  will  explore  strategies  to  improve  AP  registration  accuracy,  and 
streamline  procedures  to  enable  online  use. 
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Abstract 

Purpose:  To  adapt  a  magnetization-prepared  RF  cycled  spiral  imaging  technique,  termed  Tiprep, 
for  time-efficient  RF-insensitive  clinical  prostate  Ti  quantification  at  1.5  tesla,  and  evaluate  SNR 
limits  to  voxel-based  Ti  mapping. 

Materials  and  Methods:  A  magnetization-prepared  spiral  imaging  technique  was  adapted  for 
robust  Ti  contrast  development,  multi-slice  imaging  within  5-minutes,  and  data  regression  to  a 
monoexponential  decay.  In  vitro  testing  evaluated  RF  insensitivity  of  the  multi-slice  acquisition 
plus  method  accuracy.  A  clinical  pilot  study  was  performed  in  fifteen  patients  with  low  or 
intermediate  risk  localized  prostate  cancer. 

Results:  The  multi-slice  design  displayed  excellent  RF  insensitivity  (<1%  error  for  RF  mis- 
tunings  to  ±20%),  and  accuracy  (within  0.5%  of  gold  standard  for  Ti  values  within  1500  ms).  A 
clinical  pilot  study  reported  equivalence  between  central  gland  (CG)  and  peripheral  zone  (PZ)  Ti 
values  (CG:  1321±45  ms,  n=14;  PZ:  1359±204  ms,  n=l  1)  and  reduced  tumor  Ti  (1 196±106  ms, 
p=0.028  and  0.048  compared  to  CG  and  PZ).  SNR  measurements  identified  an 
inappropriateness  of  voxel-based  analysis. 

Conclusions:  Tiprep  can  quantify  prostate  Ti,  as  an  adjunct  measure  for  quantitative  perfusion 
measurements  in  radiotherapy  trials.  Intra-patient  heterogeneities  support  Ti  assessment  within 
individual  patients.  SNR  calculations  will  guide  protocol  improvements  towards  voxel-based 
analysis  in  future  trials. 

Keywords:  TI  quantification,  prostate  gland,  quantitative  MRI 
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Introduction 

Ti  quantification  is  motivated  primarily  by  a  need  to  quantify  exogenous  eontrast  agent 
concentration,  both  from  signal  time-eourses  and  for  molecular  targeting  experiments  (1). 

Indeed,  Ti  aecuracy  and  precision  are  acknowledged  as  prineipal  biases  for  perfusion  imaging 
analysis  (2).  Furthermore,  endogenous  Ti  retains  stand-alone  value,  for  eharaeterization  of 
partial  oxygen  pressures  plus  elucidation  of  hemorrhagic  pathology  (3,  4). 

Ti  remains  diffieult  to  measure  in  vivo,  despite  the  longevity  of  the  field.  Classieal 
inversion-recovery  approaehes  are  prohibitively  long  for  elinieal  applieation,  while  time-effieient 
and  eommereially-aecessible  methods  assume  knowledge  of  the  radio-frequeney  (RF)  field  (5, 

6).  Look-Locker  and  variable-flip-angle  standards  (VFA)  have  demonstrated  some  suceess  for 
biomedieal  translation,  but  the  RF  field  must  be  measured  independently  and  the  flip  angle 
reduetion  eompromises  signal-to-noise  ratio  (SNR)  (7). 

This  researeh  describes  the  adaption  of  a  magnetization-prepared  spiral  imaging  strategy 
for  Ti  quantification  of  the  prostate  gland,  termed  Tlprep.  The  method  provides  multi-sliee 
coverage  of  the  entire  gland  with  high  RF  insensitivity  and  a  elinically  aceeptable  time- 
efficieney  (5  minutes).  A  clinical  pilot  study  involved  15  patients  with  newly  diagnosed 
loealized  prostate  eancer  to  evaluate  zonal  and  tumor  Ti.  SNR  analysis  charaeterized  thermal 
noise  contributions  to  T i  precision,  to  guide  protocol  design  for  future  voxel-based  analysis  as  an 
adjunct  measure  for  quantitative  perfusion  measurements. 


Methods 
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All  studies  were  performed  as  part  of  an  approved  prospective  clinical  trial.  Magnetic  resonance 
(MR)  scanning  was  performed  using  a  1.5  tesla  Signa  (General  Electric  Medical  Systems, 
Waukesha,  WI). 

Methodology  development: 

Tlprep  (Fig.  1)  can  be  perceived  as  a  time-efficient  adaptation  of  a  classical  ‘single-point’ 
technique  from  T2-prepared  approach  (4),  so  that  individual  inversion  times  (TI)  are  sampled  in 
sequential  repetition  times  (TR)  with  robust  contrast  development  and  the  entire  longitudinal 
magnetization  is  sampled  using  a  90°  excitation  pulse.  To  achieve  a  clinically  realistic  scan  time 
on  the  order  of  5 -minutes,  spiral  imaging  accelerates  k- space  sampling  and  RF  cycling 
minimizes  TI  sampling  to  2  time-points  (Tf,  TE). 

Ti  contrast  is  developed  using  a  non-selective  adiabatic  inversion  pulse  followed  by  an 
arbitrary  delay  of  duration  TI.  Slice  selection  and  spatial  encoding  are  performed  using  a 
spectral-spatial  pulse  and  spiral  imaging  gradient,  with  trailing  spoiler  gradient.  Multi-slice 
imaging  is  performed  by  repetitive  slice  selection  and  spatial  encoding  at  frequent  intervals  (~ 
every  20ms  in  different  slices).  A  second  delay  of  value  TT  trails  a  final  spoiler  gradient. 
Following  TT,  a  BIR-4  adiabatic  half-passage  pulse  is  applied  to  null  the  residual  longitudinal 
and  transverse  magnetization.  By  preserving  a  constant  summation  of  TI  and  TT  across  variable 
TI,  the  period  of  longitudinal  recovery  is  kept  independent  of  TI  selection  for  any  TR. 

Central  to  Tlprep  is  RF  cycling  of  an  inversion  pulse  on/off  on  every  other  sequence 
iteration  (8,  9).  The  difference  signal  is  a  Ti -weighted  signal  decay,  while  the  additive 
component  of  Ti  relaxation  is  nulled.  TI  fitting  as  a  monoexponential  decay  rather  than  as  a 
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monoexponential  reeovery  provides  several  benefits  —  TI  sampling  ean  be  reduced  to  only  2  time 
points,  2-parameter  fitting  of  a  decay  provides  good  noise  performance,  confusion  surrounding 
magnetization  is  removed,  and  the  measurement  is  inherently  insensitive  to  imperfect  inversion 
pulses.  Because  Ti  extraction  from  Tlprep  data  depends  on  the  absolute  TI  separation  rather 
than  absolute  TI  values,  consistent  TI  values  can  be  utilized  for  all  slices  of  a  multi-slice 
acquisition. 

Quality  assurance  testing 

Demand  for  reproducible  tissue  measures  merited  quality  assurance  testing  in  static  phantoms. 
The  specific  studies  were  as  follows: 

Reproducibility  testing  The  reproducibility  of  the  multi-slice  acquisition  in  a  static  phantom  was 
tested  using  a  750ml  volume  of  water  doped  with  manganese  chloride  (MnCb)  placed  in 
longitudinal  orientation  in  a  standard  head  coil  ((Ti  ~  610  ms).  Tlprep  parameter  selection 
mapped  to  the  proposed  clinical  parameterization  (TI  of  14  and  1014  ms,  repetition  time  (TR)  of 
4400  ms,  1.1  mm  resolution  over  a  20cm  field-of-view  (FOV),  6  mm  slice  thickness,  2  nex,  5 
minute  duration).  Multi-slice  acquisitions  were  then  repeated  following  deliberate  mis-tunings 
of  the  RF  amplifier  to  almost  ±30%.  A  standard  double-angle  approach  (gradient-recalled  echo, 
flip  angles  of  60  and  120°,  TR  =  6  seconds,  128  x  128  over  20  cm  field-of-view,  single-slice) 
validated  the  RF  distribution  (10).  Measurements  were  repeated  at  all  RF  offsets  using  a 
standard  variable-flip-angle  approach  (gradient-recalled-echo  method,  TE/TR  =  4.7/1 1  ms,  0  = 
2°,  10°,  20°,  30°,  4  averages)  (5). 

Accuracy  testing  Tlprep  accuracy  was  tested  using  seven  serial  dilutions  of  Gd-DTPA 
(Magnevist,  Berlex  Canada)  in  2cc  eppendorff  tubes  (Ti  values  between  400  and  2200ms). 
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Samples  were  placed  upright  within  an  eppendorff  tube  rack  and  submerged  in  water  heavily- 
doped  with  MnC12,  to  alleviate  susceptibility  effects  at  the  air  interfaces.  Tlprep  was  applied 
with  the  following  parameters  (TI  =  14,  414,  814,  1214,  1614,  and  2014  ms;  TR  =  6000  ms;  90° 
flip  angle;  12  10.4-ms  spiral  interleaves  at  125  kHz  readout  bandwidth  over  a  16cm  field-of-view 
for  1.1mm  in-plane  resolution;  2  min  54  sec  per  Tl  value).  A  vendor-standard  inversion- 
recovery  spin-echo  sequence  provided  comparative  gold  standard  quantification  (TE=10  ms,  Tl 
=  50,  400,  800,  1200,  1600,  2000,  and  4000  ms,  128x128  matrix  over  18  cm,  50%  k-space 
acquisition,  between  10.5  and  14.5  minutes  per  TI).  TR  was  incremented  to  maintain  a  constant 
TR  and  TI  difference  of  1 1  seconds. 

Clinical  pilot  study 

A  clinical  pilot  study  involved  15  patients  with  low  or  intermediate  risk  localized  prostate  cancer 
and  no  history  of  prior  therapy.  The  patients  were  accrued  in  an  ongoing  Phase  2  clinical  trial 
designed  to  evaluate  the  impact  of  integrating  MRI  and  advanced  image-guidance  towards 
reducing  toxicity  and  improving  quality  of  life  after  high  dose  intensity  modulated  radiotherapy. 
At  least  3  days  after  ultrasound  guided  insertion  of  three-intraprostatic  gold  fiducial  markers,  all 
patients  were  scheduled  for  MRI  of  the  pelvis. 

Patients  were  positioned  in  supine  position  on  the  MRI  table  with  a  torso-phased  array 
coil  (MEDRAD  ATD  torso  device,  Warrendale,  PA)  positioned  anteriorly  and  superiorly  to  the 
prostate  gland.  Scanning  as  a  series  of  axial  acquisitions  proceeded,  including  anatomical  fast- 
spin-echo  imaging  of  the  prostate  region  (TE/TR  =  96/5000  ms,  320x256  matrix  over  20  cm), 
diffusion-weighted  echo-planar  imaging  (TE/TR  =  64/5575  ms,  128x128  matrix  over  20  cm, 
6mm  slices,  b-values  of  0  and  600  s/mm^),  and  quantitative  Ti  mapping  (TI  =  14,  1014  ms,  TR  = 
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4400  ms,  1.1x1. 1mm  resolution  over  a  20em  FOV,  6  mm  sliee  thickness,  2  nex,  5  minute 
duration).  A  stack  of  up  to  twelve  slices  provided  full  gland  coverage  (72  mm). 

Data  analysis 

Quality  Assurance  In  vitro  analysis  was  performed  using  eFilm  (Merge  Healthcare,  Milwaukee, 
WI),  with  ROIs  drawn  on  signal  images  rather  than  on  parametric  maps.  T1  fitting  to  a 
monoexponential  decay  used  Origin  software  (OriginLab  Corporation,  Northampton,  MA). 
Origin  software  was  also  utilized  for  linear  regression  in  variable-flip-angle  measurements  of  Ti. 
Testing  of  methodological  accuracy  used  standard  Bland- Altman  methodology.  Algebraic  tasks 
were  performed  using  Matlab  (The  Mathworks,  Natick,  MA). 

Prostate  Tj^  characterization  Data  analysis  involved  regions-of-interest  (ROIs)  corresponding  to 
the  central  gland  (CG),  uninvolved  peripheral  zone  (PZ),  and  tumor-dense  bearing  regions, 
delineated  as  hypointensities  on  diagnostic  fast-spin  echo  images  and  apparent  diffusion  maps  by 
an  experienced  radiation  oncologist  using  mipav  (National  Institutes  of  Health,  Bethesda,  MD, 
USA).  The  urethra  and  intra-prostatic  fiducial  markers  were  excluded  from  ROIs.  ROIs  were 
applied  to  Tlprep  signal  images  and  manually  corrected  to  account  for  evident  motion.  Ti  fitting 
used  Origin  software.  Signal-to-noise  (SNR)  calculations  used  the  methods  of  Constantinides  et 
al,  for  phased-array  coil  acquisitions  (11). 

Results: 


Quality  assurance 
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Using  a  centrally-placed  ROI  in  each  of  ten  6-mm  thick  slices,  Tlprep  preserved  the  Ti  standard 
deviation  of  the  multi-slice  acquisition  well  within  1%  for  RF  offsets  between  ±25%  (Fig.  2(a, 
b)).  Comparatively,  Ti  error  for  a  standard  VFA  approach  increased  by  approximately  20%  for 
every  10%  offset  of  the  RF  field.  Tlprep  also  demonstrated  excellent  agreement  with  the  gold- 
standard  inversion  recovery  approach  (Fig.  2(c,  d)).  Across  all  data  points,  Bland- Altman  bias 
and  limits  of  agreement  were  37  and  169  ms.  However,  the  analysis  suggested  an  increasing 
disagreement  between  measurements  as  sample  Ti  approached  2000  ms.  Ti  agreement  was 
within  27  ms  (i.e.  0.5%)  for  Ti  values  of  1500  ms  or  less. 

Clinical  pilot  study  -  Prostate  Ti  quantification 

Representative  axial  Tlprep  images  of  the  prostate  gland  at  base,  mid,  and  apical  locations  are 
presented  in  Fig.  3.  Figure  4  displays  Ti  maps  through  tumor-bearing  regions  of  the  peripheral 
zone,  plus  supporting  diagnostic  ADC  maps  and  FSE  images. 

From  clinical  acquisition,  the  mean  zonal  Ti  values  were  equivalent  at  1321±45  ms  (n  = 
14)  and  1359±204  ms  (n  =  1 1)  for  CG  and  PZ  respectively  (p  =  0.50).  Inter-patient  Ti 
heterogeneity  was  significant,  with  individual  ROIs  ranging  in  Ti  from  1237  to  1381  ms  for  CG, 
and  from  1060  to  1778  ms  for  PZ.  Ti  values  in  tumor-dense  ROIs  were  1 196±106ms  (minimum: 
1040  ms;  maximum:  1335  ms,  n  =  8).  Tumor  Ti  were  significantly  shorter  than  both  CG  Ti 
(difference  =  1 16±1 18  ms,  p=0.0276)  and  PZ  Ti  (difference  =  179±212  ms,  p=0.0484)  (Fig.  4). 

From  histogram  analysis,  95%  of  the  delineated  voxels  presented  with  SNR  values  of  32 
(CG)  and  29  (PZ)  or  more  at  TIi.  At  TU,  95%  of  SNR  values  were  greater  than  15  for  CG  and 
13  for  PZ  so  that  Ti  estimates  were  effectively  unbiased  by  the  noise  baseline  in  magnitude 


images. 


9 


Discussion 

This  research  provides  a  novel  technical  solution  for  the  challenging  problem  of  RF-insensitive 
and  time-efficient  quantitative  prostate  Ti  measurements.  The  method  leverages  the  robust 
contrast  generation  of  a  non-selective  adiabatic  inversion  pulse,  the  rapid  multi-slice  imaging 
capability  of  spiral  imaging,  and  RF  cycling  to  reduce  TI  sampling  requirements  while  providing 
an  advantageous  dependence  between  signal-to-noise  and  Ti  precision.  Clinical  investigation  in 
patients  with  low  and  intermediate  risk  prostate  cancer  demonstrated  equivalence  between  CG 
and  prostate  Ti,  and  a  significantly  reduced  tumor  Ti.  SNR  analysis  identified  a  minimum- 
useful  sampling  volume  for  thermal-noise  insensitive  Ti  quantification,  and  identified  SNR 
requirements  for  future  trials  planning  voxel-based  analysis. 

The  Tlprep  design  has  similarities  to  the  RF  pulse  train  approaches  of  Hsu  and  co¬ 
workers  because  RF  cycling  is  used  to  minimize  the  additive  contribution  to  Ti  relaxation  (8,  9). 
However,  Hsu  samples  the  resultant  Ti  decay  with  a  train  of  RF  pulses  with  reduced  RF  flip 
angle,  while  the  RF  pulse  train  in  Tlprep  is  utilized  for  multi-slice  imaging  and  a  90°  flip  angle 
is  retained  to  maximize  SNR  in  each  slice.  The  method  may  not  be  as  rapid  as  published 
applications  of  quantitative  Ti  strategies  using  variable- flip-angle  (VFA,  DESPOTl),  inversion- 
recovery  with  ss^  readout  (ir-ssfp),  and  look-locker  echo-planar  imaging  (LL-EPI)  (5,  7,  12, 

13).  However,  methodology  comparisons  are  constrained  by  the  predominant  neurology  focus  of 
these  trials,  where  the  superb  sensitivity  of  current  head  coils  can  be  leveraged  for  scan  time 
reduction.  In  addition,  variable-flip-angle  and  ss^-based  techniques  are  inherently  sensitive  to 
RF  field  heterogeneity,  which  necessitates  independent  RF  field  mapping  (5).  The  most  relevant 
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existing  comparison  may  be  the  application  by  Kershaw  et  al  of  a  3D-inversion  recovery  turbo 
field  echo  technique  for  prostate  evaluation  at  1.5  tesla,  achieving  2x2x5  mm  interpolated  voxels 
within  a  4-minute  scan  time  though  representative  images  were  not  provided  (14). 

The  high  RF  insensitivity  of  the  Tlprep  technique  should  be  a  particular  strength  of  the 
design,  and  should  have  increasing  value  as  higher  field  strength  applications  are  developed.  In 
comparison,  a  standard  VFA  approach  with  flip  angles  of  2°,  9°,  and  19°  would  introduce  Ti 
biases  of  approximately  ±20  and  ±40%  at  RF  offsets  of  ±10  and  ±20%,  while  Tlprep  maintains 
T2  biases  on  the  order  of  1%  or  less  over  a  similar  range  of  RF  mis-tuning. 

Prostate  Ti  characterization  Prostate  T 1  quantification  seems  to  be  sparsely  covered  in  the 
literature,  and  existing  publications  seem  to  have  reasonable  to  very  good  agreement  with  the 
current  results.  Noteworthy  is  Kjaer’s  early  application  of  an  inversion-recovery  approach  with 
6  inversion  times  at  1.5  T  (15).  He  identified  Ti  values  of  1022  ms  for  normal  prostate  tissue 
and  1050  ms  for  carcinoma.  Similarly,  Kershaw  et  al  measured  Ti  using  a  variable-flip-angle 
approach  and  identified  an  equivalent  normal  prostate  and  adenocarcinoma  Ti  in  2 1  patients 
(normal  prostate:  922  ms  (692  -  1130  ms);  tumor:  916  ms  (795  -  1172  ms))  (16).  More  recently, 
Kershaw  et  al.  identified  median  central  gland  and  peripheral  zone  Ti  values  of  1390  and  1360 
ms  at  1.5  tesla,  in  patients  with  benign  prostatic  hyperplasia  and  an  absence  of  prostate  cancer. 
The  intra-quartile  ranges  were  approximately  10%  (14).  In  a  prior  study,  DeBazelaire  et  al. 
measured  mean  and  standard  deviation  prostate  gland  Ti  values  of  1317±85  ms  at  1.5  tesla  and 
1597±42  ms  at  3.0  tesla  using  matched  regions-of-interest  (17).  Noworolski  et  al.  quantified 
prostate  Ti  values  at  1.5  T  using  a  spin-echo  method  at  up  to  four  repetition  times,  and  reported 
Ti  of  15701622  ms  for  normal  peripheral  zone  and  16701421  ms  for  peripheral  zone  cancer  (18). 
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SNR  analysis  SNR  estimates  are  signifieant  because  they  define  the  best  case  precision  of  the  Ti 
measurement.  Based  on  Monte  Carlo  simulation,  a  monoexponential  decay  at  experimental  Tl 
and  Tl  of  1300  ms  identifies  a  noise-mediated  Ti  standard  deviation  (aTi)  of  5%,  2.5%,  2%, 
1.5%,  and  1%  at  SNR  values  of  35,  70,  90,  115,  and  170.  oTi  for  fitting  a  2-parameter 
monoexponential  recovery  with  3  TI  and  equivalent  Tl  and  SNR  is  similar,  but  care  must  be 
taken  to  judge  the  polarity  of  small  magnitude  signals. 

Using  Tlprep,  histogram  analysis  reported  a  SNR  of  32  or  more  in  at  least  95%  of  voxels 
at  Til,  so  that  voxel-based  analysis  would  be  prone  to  variability  from  thermal  noise  (oTi  ~  5% 
at  SNR  of  35).  One  may  then  combine  voxels  into  a  ROI  to  increase  robustness  to  thermal  noise, 
and  the  SNR  gain  is  proportional  to  the  square  root  of  the  number  of  independent  voxels  within 
the  ROI.  For  example,  combination  of  10  voxels  (i.e.  73  mm  sampling  volume)  with  SNR  30 
would  maintain  thermal  noise  variability  within  2%,  while  a  35  voxel  ROI  (254  mm  sampling 
volume)  would  maintain  thermal  noise  variability  within  1%.  Consequently,  noise  precision  is 
traded  off  for  volume  averaging  in  the  absence  of  other  strategies  to  augment  SNR. 

It  is  worthwhile  to  note  that  TI  selection  can  be  tailored  for  the  target  Ti.  Roughly,  to 
maximize  precision  of  monoexponential  fitting  using  two  data-points  only,  Tf  should  be 
minimized  while  TU  should  be  1.1  times  the  expected  time-constant  (19).  For  this  pilot  trial,  TU 
was  1014  ms  so  that  an  incremental  increase  in  noise  performance  could  be  gained  via 
prolongation  of  TU  towards  1400  ms  at  experimentally  measured  Ti.  Slice  number  may  need  to 
be  constrained  for  contrast-enhanced  studies  because  T i  decay  during  the  imaging  interval  may 
limit  SNR  independently  of  TI  selection. 
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Clinical  significance  The  primary  clinical  significance  of  T i  quantification  is  its  central  role  as 
an  adjunet  measure  for  quantitative  analysis  of  dynamic  contrast  enhancement  (DCE)  (2,  20). 
These  results  highlight  a  eonsiderable  intra-patient  zonal  and  T  i  heterogeneity  whieh  further 
supports  individual  Ti  assessment  in  DCE  examinations.  Ideally,  a  progression  of  dynamic 
prostate  evaluation  to  voxel-based  analysis  will  be  eoncomitant  with  a  progression  of 
quantitative  Ti  evaluation  within  registered  voxels.  Part  and  pareel  of  this  progression  to  voxel- 
based  parametric  mapping  will  be  a  growing  importanee  of  thermal  noise  variability,  particularly 
for  regression  analysis  with  sparse  samples  (i.e.  ADC,  Ti,  T2).  Practical  solutions  for 
measurements  with  noise-limited  precision  at  constant  voxel  volume  and  scan  time  include 
improved  RE  eonfigurations  (i.e.  endo-rectal  eoils  in  tandem  with  torso  phased-arrays)  and 
stronger  magnets. 

Conclusions 

A  magnetization-prepared  spiral  imaging  strategy  has  been  adapted  for  time  efficient  and  RF- 
insensitive  prostate  T 1  quantification.  Applied  to  patients  with  low  or  intermediate  risk  localized 
prostate  eancer  and  no  treatment  history,  it  demonstrated  equivalence  between  mean  zonal  Ti 
values  and  a  reduced  tumor  Ti.  However,  intra-patient  zonal  and  tumor  Ti  heterogeneities  were 
considerable  so  that  individual  prostate  Ti  mapping  is  merited  for  quantitative  assessment  of 
contrast  enhancement  dynamics.  SNR  analysis  identified  eonstraints  for  thermal  noise 
insensitive  Ti  measurements,  to  guide  protoeol  design  for  voxel-based  Ti  quantification  as  an 
adjunct  measure  to  support  voxel-based  DCE  analysis. 
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Figure  Legends: 

Figure  1:  Tlprep  pulse  sequence  design.  Ti  contrast  is  developed  globally  during  a  preparation 
interval  of  duration  Tl,  termed  the  inversion  time.  Multi-slice  imaging  is  performed  using  a 
spectral- spatial  pulse  and  spiral  imaging  gradients,  applied  repetitively  for  each  imaging  slice 
(‘N’  times  for  ‘N’  slices).  Following  a  duration  TF,  an  adiabatic  BlR-4  pulse  is  applied  to  null 
the  remnant  longitudinal  magnetization.  Measurement  independence  on  the  TR  is  achieved  by 
maintaining  a  consistent  summed  Tl  -I-  TF  for  all  inversion  times.  On  every  other  sequence 
iteration,  RF  cycling  of  a  second  adiabatic  inversion  pulse  is  applied  prior  to  the  start  of  the 
imaging  interval  to  null  the  additive  Ti  recovery  term. 

Figure  2:  Quality  Assurance  testing:  (a)  Reproducibility  phantom,  displaying  B1  map  along  a 
coronal  section  of  a  MnC12-doped  water  bottle  phantom  (long  rectangle),  and  location  of  10 
imaging  slices  plus  approximate  diameter  of  ROIs  in  Tlprep  images  (short  parsed  rectangle);  (b) 
Percentage  Ti  error  as  a  function  of  RF  offset,  demonstrating  negligible  RF  sensitivity  using 
multi-slice  Tlprep  and  considerable  RF  sensitivity  using  a  variable-flip-angle  approach;  (c) 
Tlprep  images  of  an  accuracy  phantom,  consisting  of  7  eppendorff  tubes  doped  to  variable 
extent  with  Gd-DTPA,  and  oriented  vertically  in  a  bath  of  water  heavily  doped  with  MnCb;  (d) 
Bland-Altman  analysis  demonstrating  equivalence  between  Tlprep  and  the  gold  standard 
inversion-recovery  technique.  The  bias  and  limit  of  agreement  (LOA)  are  displayed. 

Figure  3:  Representative  clinical  images  at  three  levels  of  the  prostate  gland:  (top  row;  a-c) 
standard  diagnostic  FSE  sections;  (middle  row;  d-f)  Tli  Tlprep  images;  (bottom  row;  g-i)  TF 
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Tlprep  images.  Implanted  gold  intra-prostatic  fiducial  marker  are  visible  as  focal 
hypointensities  in  in  each  spiral  image. 

Figure  4  -  Representative  images  from  three  patients,  termed  patients  A,  B.  and  C;  (left)  ADC 
map;  (mid)  diagnostic  FSE  image;  and  (right)  quantitative  Ti  maps  (windowed  between  0  and 
2000ms).  For  each  cases,  hypo-intense  regions  deemed  to  be  tumor-dense  can  be  visualized  in 
both  ADC  and  FSE  acquisitions.  Patient  A  displayed  Ti  values  of  1272  ms  (CG),  1349  ms  for 
peripheral  zone,  and  1 162  ms  for  tumor.  Patient  B  displayed  Ti  values  of  1329  ms  (CG),  1311 
ms  for  peripheral  zone,  and  1 174  ms  for  tumor. 
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Figure  1:  Tlprep  pulse  sequence  design.  Ti  contrast  is  developed  globally  during  a  preparation 
interval  of  duration  TI,  termed  the  inversion  time.  Multi-slice  imaging  is  performed  using  a 
spectral-spatial  pulse  and  spiral  imaging  gradients,  applied  repetitively  for  each  imaging  slice 
(‘N’  times  for  ‘N’  slices).  Following  a  duration  TF,  an  adiabatic  BIR-4  pulse  is  applied  to  null 
the  remnant  longitudinal  magnetization.  Measurement  independence  on  the  TR  is  achieved  by 
maintaining  a  consistent  summed  TI  +  TF  for  all  inversion  times.  On  every  other  sequence 
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iteration,  RF  cycling  of  a  second  adiabatic  inversion  pulse  is  applied  prior  to  the  start  of  the 
imaging  interval  to  null  the  additive  Ti  recovery  term. 
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Figure  2:  Quality  Assurance  testing:  (a)  Reproducibility  phantom,  displaying  B1  map  along  a 
coronal  section  of  a  MnC12-doped  water  bottle  phantom  (long  rectangle),  and  location  of  10 
imaging  slices  plus  approximate  diameter  of  ROIs  in  Tlprep  images  (short  parsed  rectangle);  (b) 
Percentage  Ti  error  as  a  function  of  RF  offset,  demonstrating  negligible  RF  sensitivity  using 
multi-slice  Tlprep  and  considerable  RF  sensitivity  using  a  variable-flip-angle  approach;  (c) 
Tlprep  images  of  an  accuracy  phantom,  consisting  of  7  eppendorff  tubes  doped  to  variable 
extent  with  Gd-DTPA,  and  oriented  vertically  in  a  bath  of  water  heavily  doped  with  MnCb;  (d) 
Bland- Altman  analysis  demonstrating  equivalence  between  Tlprep  and  the  gold  standard 
inversion-recovery  technique.  The  bias  and  limit  of  agreement  (LOA)  are  displayed. 
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Figure  3:  Representative  elinical  images  at  three  levels  of  the  prostate  gland:  (top  row;  a-e) 
standard  diagnostic  FSE  sections;  (middle  row;  d-f)  Tli  Tlprep  images;  (bottom  row;  g-i)  TI2 
Tlprep  images.  Implanted  gold  intra-prostatic  fiducial  marker  are  visible  as  focal 
hypointensities  in  in  each  spiral  image. 


22 
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Figure  4  -  Representative  images  from  three  patients,  termed  patients  A,  B.  and  C;  (left)  ADC 
map;  (mid)  diagnostic  FSE  image;  and  (right)  quantitative  Ti  maps  (windowed  between  0  and 
2000ms).  For  each  cases,  hypo-intense  regions  deemed  to  be  tumor-dense  can  be  visualized  in 
both  ADC  and  FSE  acquisitions.  Patient  A  displayed  Ti  values  of  1272  ms  (CG),  1349  ms  for 
peripheral  zone,  and  1 162  ms  for  tumor.  Patient  B  displayed  Ti  values  of  1329  ms  (CG),  1311 


ms  for  peripheral  zone,  and  1 174  ms  for  tumor. 
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Abstract 


T2  quantification  may  augment  diagnostic  T2-weighted  imaging;  to  improve  cancer  detection  via 
auto-segmentation  of  multi-parametric  acquisitions  and  to  potentiate  longitudinal  studies  of 
prostate  cancer.  However,  robust  quantitative  techniques  are  not  generally  accessible  or  refined 
for  clinical  translation.  This  research  describes  the  adaptation  of  a  magnetization-prepared  spiral 
imaging  technique,  termed  T2prep  for  prostate  application,  providing  whole  gland  coverage 
within  a  5-minute  interval  with  considerable  insensitivity  to  radio-frequency  (RF) 
inhomogeneities.  Clinical  piloting  in  two  cohorts  with  distinct  clinical  histories  demonstrated 
the  anticipated  differences  in  zonal  and  tumor  T2,  including  tumor  T2  shortening  compared  to 
peripheral  zone,  and  post-radiotherapy  shortening  of  peripheral  zone  T2.  SNR  calculations  were 
performed  for  data  acquired  with  or  without  an  endo-rectal  coil  in  tandem  with  a  torso  phased 
array,  to  judge  the  potential  for  voxel-based  T2  mapping  and  thereby  support  focal  biological 
characterization  of  cancer,  hypoxia,  and  response  to  therapy  within  regions  of  dense  cancer 
burden. 
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Introduction 


Quantitative  MRI  biomarkers  (i.e.  apparent  diffusion  coefficient,  vascular  metrics,  metabolite 
profiles)  are  particularly  relevant  for  cancer  management,  both  for  differential  diagnosis  and  to 
judge  the  efficacies  of  novel  therapeutic  strategies  over  time  in  pre-clinical  and  clinical  trials  (1). 
T2  relaxation  is  also  central  to  MRI  diagnostics,  and  T2  quantification  may  complement  both 
diagnostics  and  trials  of  tumor  response  (2-4).  However,  T2  quantification  is  often  restricted  to 
research  purposes  despite  the  long  and  rich  history  of  MR  relaxometry. 

Robust  T2  quantification  in  vivo  constrains  pulse  sequence  design  to  limit  biases  from 
system  imperfections  and  mis-tunings,  physiologic  motions,  and  volume  averaging,  without 
compromising  time-efficiency  and  volumetric  coverage  targets  (5,  6).  Within  this  paradigm,  a 
magnetization-prepared  spiral  imaging  technique,  termed  T2prep,  was  refined  to  provide  an 
imaging  solution  for  vascular  oximetry  measurements  and  peripheral  venous  angiography, 
though  characteristics  of  robust  global  contrast  development  and  longitudinal  storage  of  contrast 
also  support  quantitative  tissue  evaluation  (7,  8). 

This  paper  describes  an  adaption  of  T2prep  for  clinical  prostate  measurements  to  quantify 
zonal  and  tumor  T2  and  SNR  values.  Clinical  piloting  used  two  patient  cohorts,  which  are 
distinct  based  on  patient  history  of  radiotherapy  and  on  inclusion  of  an  endo-rectal  coil  for  MR 
signal  reception.  SNR  analysis  estimated  noise  precision  in  voxel-based  maps,  to  guide  protocol 
design  for  future  evaluation  of  hypoxia  and  therapy  response. 

Methods 
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All  patients  provided  informed  consent  and  were  prospectively  enrolled  in  approved  clinical 
trials.  All  studies  were  performed  using  a  1.5  tesla  MRI  scanner  (Signa,  General  Electric 
Medical  Systems,  Waukesha,  WI). 

Patient  selection 

Cohort  A 

Cohort  A  included  16  men  with  low/intermediate  risk  localized  prostate  cancer  and  accrued  in  an 
ongoing  Phase  2  clinical  trial.  At  least  3  days  after  ultrasound  guided  insertion  of  three  intra¬ 
prostatic  gold  fiducial  markers,  all  patients  were  scheduled  for  radiation  planning  CT  scanning 
(2mm  axial  cuts,  no  spacing)  followed  by  MRI  of  the  pelvis. 

Cohort  B 

Cohort  B  included  8  men  suspected  to  have  local  recurrence  following  external  beam 
radiotherapy  and  accrued  in  an  ongoing  clinical  trial  of  MRI  guided  biopsy  for  evaluation  of 
local  recurrence. 

T2Prep  acquisition  method 

Figure  1  presents  the  template  for  the  T2prep  technique  (5-8).  T2  contrast  is  prepared  using  a 
non-selective  robust  refocusing  train,  returned  to  the  longitudinal  axis.  Slice  selection  and 
spatial  encoding  are  performed  using  a  spectral-spatial  pulse  and  spiral  imaging  gradient,  with 
trailing  spoiler  gradient.  Multi-slice  imaging  is  performed  by  repetitive  slice  selection  and 
spatial  encoding  at  frequent  intervals  (~  every  20ms  in  different  slices).  A  BIR-4  adiabatic  half¬ 
passage  pulse  trails  the  final  imaging  interval  to  null  the  residual  longitudinal  and  transverse 
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magnetization.  The  delay  between  the  spiral  readout  and  BIR-4  pulse  is  shifted  to  keep  the 
period  of  longitudinal  recovery  independent  of  TE  for  the  interleaved  acquisition.  RF  cycling  is 
applied,  consisting  of  an  adiabatic  inversion  pulse  prior  to  the  first  imaging  interval  of  every 
other  sequence  iteration.  Subtraction  isolates  the  T2-weighted  decay  function  without  bias  from 
additive  Ti  relaxation  during  the  period  of  longitudinal  storage. 

Experiment!:  Quality  Assurance 

Quality  testing  evaluated  the  RF  sensitivity  of  the  multi-slice  T2prep  acquisition,  using  a  750  ml 
bottle  of  MnC12-doped  water  (T2~72  ms),  the  head  coil  of  the  GE  Signa,  and  a  12-slice 
prescription  covering  7.2  cm  in  6-mm  sections.  The  T2prep  parameter  set  included  2  echo  times 
(TE)  of  3.2  and  88.4  ms,  4400  ms  repetition  time  (TR),  l.lxl.lx6-mm  voxels  based  on  spiral  k- 
space  sampling  at  125  kHz  over  a  20-cm  field-of-view,  and  2  averages  for  a  5-minute  scan 
duration.  RF  tuning  was  validated  using  the  double-angle  technique  (gradient-recalled  echo 
acquisitions  at  flip  angles  of  60  and  120°  with  minimized  echo  times  of  1.3  ms,  long  repetition 
times  of  8500  ms,  and  128x64  matrix  over  a  20  cm  field-of-view)  (9).  The  Bi  profile  was 
calculated  after  averaging  10  adjacent  voxels  at  each  point  along  the  length  of  the  phantom.  The 
T2prep  acquisition  was  repeated  at  5  unit  increments  of  the  RF  amplifier  (ranging  from  ±20  %) 
via  deliberate  system  mis-tunings,  to  encompass  the  characteristic  ±15  %  offsets  of  a  1.5  tesla 
whole-body  scanner  (7). 

Experiment  2:  T2  quantification  in  patients  with  low  or  intermediate  risk  localized  prostate 
cancer  prior  to  external  beam  radiotherapy 
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Cohort  A  patients  were  positioned  supine  on  the  MRI  table  with  a  torso-phased  array  coil 
(MEDRAD  ATD  torso  device,  Warrendale,  PA)  placed  anterior  and  posterior  to  the  pelvis.  The 
MRI  examination  proceeded,  as  a  series  of  axial  scans  of  the  prostate  with  fixed  graphical 
prescription  and  common  20  cm  field-of-view,  including  2D-FSE  imaging  (TE/TE  =  96/5000 
ms,  echo-train  length  =  16,  320x256  matrix,  3  mm  slice  thickness),  diffusion-weighted  echo- 
planar  imaging  (TE/TR=64/5575  ms;  b  =  0  and  600  s/mm^,  128x128  matrix,  6  mm  slice 
thickness),  and  T2prep  (TE  =  3.2,  88.4  ms;  TR=  4400  ms,  1. 1x1.1  mm  voxels,  6  mm  slice 
thickness,  5  minute  scan).  The  majority  of  acquisitions  required  10  to  12  6-mm  slices  for  full 
prostate  coverage. 

Experiment  3:  T2  quantification  in  patients  with  suspected  local  recurrence  after  external 
beam  radiotherapy 

Cohort  B  patients  were  positioned  supine  on  the  MRI  table  and  an  endorectal  coil  (MEDRAD 
MRInnervu  endorectal  coil)  -  modified  to  attach  to  a  locking  device,  was  inserted.  A  pelvic 
phased  array  coil  (MEDRAD  ATD  torso  device)  was  positioned  anterior  and  posterior  to  the 
pelvis.  T2prep  (TE  =  3.2,  88.4  ms;  TR=  4400  ms,  1.1x1. 1x6  mm  voxels  over  a  20  cm  field-of- 
view,  5  minute  scan,  10-12  slices  for  typical  whole  gland  coverage)  was  then  incorporated  into  a 
multi-parametric  MRI  protocol,  which  included  2D-FSE  and  DW-EPI.  The  latter  acquisitions 
shared  similar  parameters  with  the  design  of  Experiment  2. 

Data  analysis 

For  experiment  1,  the  RF  field  was  mapped  via  standard  analysis  of  the  double-angle  method 
(Matlab,  The  Mathworks,  Natick,  MA)  (9).  Echo  time  and  signal  magnitude  data  pairs  were 


6 


extracted  using  contours  drawn  on  eFilm  (Merge  Healthcare,  Milwaukee,  WI).  T2  regression 
analyses  were  performed  using  the  Origin  software  package  (OriginLab,  Northampton,  MA). 

For  experiments  2  and  3,  signal  magnitudes  were  extracted  from  regions-of-interest 
(ROIs)  representing  the  central  gland  (CG)  and  peripheral  zone  (PZ),  which  were  delineated 
manually  using  MIPAV®  (NIH  version  4.2.1)  on  FSE  images  by  an  experienced  radiation 
oncologist.  Tumor  ROIs  were  demarcated  on  FSE  images,  guided  by  visualization  of 
characteristic  nodular  hypointensities  in  diffusion-weighted  and/or  FSE  images,  rapid  and 
intense  contrast  enhancement  and  washout  on  DCE,  and  needle  core  biopsy  results.  FSE  ROIs 
were  copied  onto  T2prep  images,  and  manually  adjusted  to  account  for  motion-related  mis¬ 
registration  and  to  exclude  artifacts  of  implanted  gold  seeds.  Signal-to-noise  for  a  given  ROI 
was  calculated  as  per  methods  of  Constantinides  et  al,  for  phased-array  coil  acquisitions  (10). 
All  significance  testing  used  the  Student’s  two-tailed  t-test,  with  significance  defined  by 
reduction  of  the  p-value  below  0.05. 

Results 

Experiment  1:  Quality  Assurance 

Figure  2(a)  displays  the  RF  uniformity  along  the  MnC12-doped  water  phantom,  as  demarcated  by 
the  solid  line.  The  dashed  box  demarcates  approximately  the  prescription  of  the  12-slice 
acquisition  plus  centrally-placed  ROIs.  Within  this  volume,  the  RF  amplitude  remained  within 
0.2  %  of  the  target  amplitude.  Figure  2(b)  displays  the  mean  and  standard  deviation  of  T2  across 
all  slices  at  a  given  RF  offset,  using  a  centrally  placed  ROI  of  3cm^.  The  T2  offset  remains  well 
within  2  %  for  RF  offsets  of  10  %  or  less,  but  elevates  towards  3  %  for  RF  offsets  of  15%.  T2 
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error  accrues  more  rapidly  at  RF  offsets  of  1 5  %  or  greater,  consistent  with  the  limits  for  success 
of  the  MLEV  pattern  of  90xl80y90x  composite  refocusing  pulses  (6). 

Experiment  2:  T 2  quantification  in  patients  with  low  or  intermediate  risk  localized  prostate 
prior  to  external  beam  radiotherapy 

Figure  3  displays  representative  T2prep  images  and  corresponding  T2  maps  across  base,  mid,  and 
apical  regions  of  the  prostate.  The  locations  of  gold  intra-prostatic  fiducial  markers  are 
displayed  as  focal  hypointensities  in  the  proton-density-weighted  images. 

Zonal  T2  values  were  78±7  ms  for  CG  and  109±25  ms  for  PZ  (p  =  0.0002).  Peripheral 
zone  tumors  were  identified  in  1 1  of  16  subjects.  At  a  mean  of  81±10  ms,  tumor  T2  was  reduced 
significantly  from  uninvolved  peripheral  zone  by  30±24  ms  (range:  9  to  91  ms,  p  =  0.0026). 
From  histogram  analysis,  95%  of  the  delineated  voxels  in  zonal  ROIs  presented  with  SNR  values 
of  28  (CG)  and  27  (PZ)  or  more,  at  the  first  TE  of  3  ms.  The  corresponding  lower  confidence 
intervals  for  SNR  at  the  second  TE  of  88ms  were  approximately  9  for  both  CG  and  PZ. 

Experiment  3:  T2  quantification  in  patients  with  suspected  local  recurrence  after  external 
beam  radiotherapy 

Representative  T2-weighted  images  and  T2  maps  at  basal,  middle,  and  apical  regions  of  the 
prostate  gland  are  illustrated  in  Fig.  4.  A  comparison  of  Cohort  A  and  B  images  in  tumor¬ 
bearing  regions  are  presented  in  Fig.  5.  Artifacts  of  the  patient  preparation  included  a  bright 
intra-urethral  signal  due  to  foley  catheterization  (3/8  patients)  and  occasional  intra-prostatic 
edema  from  injection  of  local  anesthetic  (5/9  patients).  These  regions  were  excluded  from  ROIs. 
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The  mean  T2  values  were  72±6  ms  for  the  CG,  and  94±12ms  for  the  PZ  (p  =  0.003,  n  = 
8).  Tumor-dense  regions  were  identified  in  5  of  8  subjects  (Fig.  5(b)).  The  mean  tumor  T2  was 
81±4  ms,  which  approached  significance  in  difference  with  the  uninvolved  PZ  (p=0.08,  n=5). 
Consistent  with  the  sensitivity  profile  of  the  endo-rectal  coil,  histogram  analysis  identified  a  non¬ 
normal  distribution  of  image  SNR,  so  that  voxels  had  mean  and  standard  deviation  SNRs  of 
96±21  for  CG  and  147±32  for  PZ,  but  95%  of  the  CG  voxels  had  SNR  greater  than  62  and  95% 
of  PZ  voxels  had  SNR  greater  than  98. 

Discussion 


This  research  describes  the  adaptation  of  a  T2-prepared  spiral  imaging  method  for  T2 
quantification  in  prostate  radiotherapy.  The  method  effectively  combines  RF  insensitivity  with 
time-efficiency  and  multi-slice  coverage.  Clinical  piloting  used  two  cohorts  at  distinct  stages  of 
cancer  progression  and  radiation  treatment.  Our  results  are  consistent  with  the  known  variations 
in  prostate  T2,  including  elevated  values  in  PZ  compared  to  CG  and  tumor,  and  reduced  prostate 
T2  contrast  following  external  beam  irradiation.  Specifically,  Cohort  B  presented  with  lower  PZ 
T2  at  near-equivalent  CG  and  tumor  T2  compared  to  Cohort  A.  SNR  analysis  characterized  the 
potential  for  voxel-based  T2  mapping  when  a  torso  phased  array  coil  is  used  with  and  without  an 
endo-rectal  coil.  This  standard  is  a  step  in  the  progression  towards  robust  voxel-based  T2 
mapping 

Technical  relevance 
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Motivation  to  review  T2  measurement  development  is  a  surprisingly  broad  divergence  in  prostate 
measurements.  Divergent  results  may  stem  from  the  lack  of  a  standardized  technique,  including 
factors  affecting  robustness  and  echo  time  sampling,  given  accumulating  evidence  of  non¬ 
monoexponential  prostate  T2  relaxation  (11,  12).  For  example,  Gibbs  et  al  demonstrated  zonal 
T2  of  88±13  ms  (CG)  and  122±34  ms  (PZ)  and  sampled  to  1 12  ms  (2),  Danger  et  al  demonstrated 
PZ  and  tumor  T2  of  108  and  82  ms  and  sampled  to  90  ms  (4),  Kershaw  et  al  demonstrated 
median  CG  and  PZ  T2  on  the  order  of  150  ms  and  sampled  to  240  ms  (13),  while  Roebuck  et  al 
demonstrated  suspected  healthy  and  suspected  cancer  T2  values  of  193±49  ms  and  100±26  ms 
and  sampled  to  180  ms  (14).  Current  results  are  consistent  in  TE  selection  and  T2  to  Danger  et  al 

(4). 

In  design,  T2prep  presents  with  advantages  over  the  standard  Carr-Purcell-Meiboom-Gill 
(CPMG)  refocusing  train,  because  it  uses  a  validated  combination  of  non-selective  refocusing 
pulses  and  scan  time  does  not  limit  slice  coverage.  However,  scan  time  increases  linearly  with 
the  number  of  echo  times,  which  constrains  T2  fitting  to  a  monoexponential  decay  if  only  2  echo 
times  are  acquired.  This  paper  presented  a  5  minute  whole-gland  acquisition  with  7.3  mm^ 
voxels.  In  comparison,  Roebuck  et  al.  demonstrated  a  10.7-minute  voxel-based  T2  mapping  of 
the  prostate  gland  with  0.55  x  1.09  x  4  mm  voxels  (2.4  mm^  volume)  using  a  CPMG  technique 
augmented  by  an  endo-rectal  coil  in  combination  with  a  4-coil  torso  phased  array  at  1.5  tesla 
(14).  Kershaw  et  al  utilized  parallel  imaging  to  reduce  the  data  acquisition  time  of  a  CPMG 
approach  to  4-minutes,  using  a  cardiac  phased-array  coil  at  1.5  tesla  with  2.3  mm  interpolated  in¬ 
plane  resolution  and  a  4.5  mm  slice  thickness  (at  least  23.8  mm^  voxels)  (13).  However,  CPMG 
refocusing  pulses  must  be  slice-selective  to  gain  the  time  efficiency  of  slice-interleaving  for 
multi-slice  acquisitions,  which  artificially  accelerates  T2  relaxation  (15). 
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SNR 


By  definition,  SNR  defines  the  maximal  precision  for  T2  quantification.  SNR  analysis  is 
therefore  requisite  for  progression  from  ROI  to  voxel-based  T2  mapping,  to  support  diagnostics 
based  on  auto-segmentation,  biopsy  guidance,  tumor  hypoxia  quantification,  and  tumor  radiation 
response  quantification. 

From  Monte  Carlo  simulation  of  a  monoexponential  decay  at  experimental  TE  and  T2  of 
80ms,  T2  standard  deviation  from  thermal  noise  reduces  below  5%,  2.5%,  2%,  1.5%,  and  1%  for 
SNR  values  of  35,  70,  90,  115,  and  170.  From  histogram  analysis  of  Cohort  A,  95%  of  the 
voxels  presented  with  SNR  values  of  28  (CG)  and  27  (PZ)  or  more,  at  the  first  TE  of  3  ms. 
Consequently,  any  voxel-based  analysis  would  be  quite  prone  to  thermal  noise  variability.  For 
cohort  B,  95%  of  the  voxels  in  zonal  ROIs  presented  with  SNR  values  of  62  (CG)  and  98  (PZ). 
Consequently,  voxel-based  analysis  of  the  T2prep  method  approaches  feasibility  when  an  endo¬ 
rectal  coil  is  used  in  tandem  with  a  torso-phased  array. 

At  all  TE  values,  SNR  must  also  be  significantly  greater  than  the  noise  magnitude,  else 
the  T2  estimate  is  artificially  elevated  by  the  noise  baseline  signal.  For  Cohort  A,  95%  of  voxels 
had  SNR  of  greater  than  9  at  the  second  TE  of  88ms,  so  that  noise  bias  was  not  significant . 

Clinical  relevance 

T2-weighted  imaging  provides  moderate  sensitivity  (50-77%)  and  specificity  (36-63%)  for 
identifying  tumor-dense  regions  of  the  prostate  gland,  with  improved  diagnostic  performance  for 
combined  pelvic  phased  array  and  endo-rectal  coil  acquisitions  (77-96%)  (16-17).  However, 
diagnostic  performance  is  reduced  in  patients  with  small  tumor  burden  (size<  10  mm,  volume<5 
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ml)  (17),  transition  zone  or  anterior  disease  (16),  low  tumor  density  (5),  or  prior  history  of 
external  beam  irradiation  (18).  Diagnostic  performance  is  improved  when  T2- weighted  images 
are  reviewed  in  parallel  with  images  weighted  by  water  diffusion,  dynamic  contrast 
enhancement,  and/or  *H-spectroscopy  (1-4,19-20).  T2  quantification  may  further  augment  a 
multi-parametric  analysis  of  prostate  cancer,  by  potentiating  auto-segmentation  of  tumor-dense 
regions  on  parametric  maps  (T2,  DCE,  ADC)  and  for  more  comprehensive  longitudinal 
monitoring  of  tumor  response  to  therapy. 

Conclusion 

A  magnetization-prepared  spiral  imaging  technique  was  adapted  for  prostate  T2  quantification 
with  adequate  time  efficiency,  volumetric  coverage,  and  RF  insensitivity  for  clinical  piloting. 
Consistent  with  prior  literature,  cohorts  without  history  of  external  beam  irradiation 
demonstrated  significantly  prolonged  T2  relaxation  in  the  uninvolved  peripheral  zone  compared 
to  the  central  gland  zone  and  tumor.  Regional  T2  differences  were  reduced  following  external 
beam  radiation.  SNR  must  be  augmented  from  that  provided  by  standard  phased-array  coil 
acquisitions  to  guide  designs  of  future  trials  targeting  focal  tumor  hypoxia  and  therapy  response 
based  on  voxel-based  T2  mapping. 
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Figure  Legend: 


Figure  1 :  T2prep  design,  consisting  of  sequential  intervals  for  contrast  preparation,  imaging,  and 
nulling  of  remnant  magnetization.  Contrast  preparation  uses  composite  90  and  -90  pulses, 
bracketed  around  a  MLEV  train  of  composite  90xl80y90x  refocusing  pulses  (4  shown),  and  a 
trailing  spoiler  pulse  and  gradient.  An  adiabatic  inversion  pulse  is  applied  on  every  even- 
numbered  excitation  to  preserve  T2  contrast  during  longitudinal  storage.  Imaging  uses  a 
spectral-spatial  excitation  and  spiral  imaging  gradients,  with  trailing  spoiler  gradient.  A  BIR-4 
pulse  and  trailing  spoiler  resets  the  magnetization  to  ensure  constant  period  of  T1  recovery  for  all 
echo  times. 

Figure  2:  Quality  assurance  testing;  (a)  Coronal  B1  map  through  the  phantom.  The  height  of  the 
narrow  rectangle  demarcates  the  volume  averaged  for  B1  quantification,  while  the  dashed  box 
demarcates  the  location  of  T2  evaluation.  The  length  of  the  dashed  box  demarcates  the  slice 
package,  while  the  height  of  the  dashed  box  demarcates  the  diameter  of  the  ROI  for  noise- 
insensitive  T2  quantification,  (b)  T2  variation  across  sequential  imaging  slices  at  RF  field  offset 
between  approximately  -20  and  20%.  The  T2  error  is  the  percentage  standard  deviation  of  T2 
measurements  across  all  slices,  relative  to  the  T2  value  in  the  central  slice  at  minimal  RF  offset. 

Figure  3:  T2prep  acquisitions  and  comparative  FSE  images  acquired  from  a  Cohort  A  subject 
throughout  the  prostate;  column  (A)  Early  echo  time  T2prep  image;  column  (B)  late  echo  time 
images;  and  column  (C)  FSE  image.  The  T2prep  images  are  windowed  between  0  and  1200  and 
0  and  500  signal  units  respectively.  Similarly,  FSE  images  are  windowed  between  0  and  500 


16 


signal  units.  Intra-prostatic  implanted  gold  fiducial  markers  are  visible  as  focal  hypo-intensities 
in  some  figure  panels  (white  arrows  in  column  (A)). 

Figure  4:  Representative  T2prep  acquisitions,  acquired  at  apex,  mid  (x2),  and  base  locations  of 
the  prostate  using  an  endo-rectal  coil;  Column  (A)  Early  echo  time  images,  windowed  between 
signal  intensities  of  0  and  1000  (256  grey  levels);  Column  (B)  late  echo  time  images,  windowed 
between  signal  intensities  of  0  and  500  (256  grey  levels);  and  (c)  T2  maps,  windowed  between 
T2  values  of  0  and  200ms  (256  grey  levels).  Three  intra-prostatic  implanted  gold  fiducial 
markers  are  visible  as  focal  hypo-intensities  in  Column  (A),  and  demarcated  with  the  black 
arrows. 

Figure  5:  Representative  images  for  Cohort  A  (top  rows)  and  Cohort  B  (bottom  rows), 
displaying  tumors  as  nodular  hypointensities  in  diagnostic  fast-spin-echo,  T2-weighted  T2prep 
images,  and  T2  maps  created  from  T2prep  images  at  echo  times  of  3  and  88ms.  The  T2  maps 
are  windowed  between  0  and  200ms. 
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Figure  1:  T2prep  design,  consisting  of  sequential  intervals  for  contrast  preparation,  imaging,  and 
nulling  of  remnant  magnetization.  Contrast  preparation  uses  composite  90  and  -90  pulses, 
bracketed  around  a  MLEV  train  of  composite  90xl80y90x  refocusing  pulses  (4  shown),  and  a 
trailing  spoiler  pulse  and  gradient.  An  adiabatic  inversion  pulse  is  applied  on  every  even- 
numbered  excitation  to  preserve  T2  contrast  during  longitudinal  storage.  Imaging  uses  a 
spectral-spatial  excitation  and  spiral  imaging  gradients,  with  trailing  spoiler  gradient.  A  BIR-4 
pulse  and  trailing  spoiler  resets  the  magnetization  to  ensure  constant  period  of  T1  recovery  for  all 
echo  times. 
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Figure  2:  Quality  assurance  testing;  (a)  Coronal  B1  map  through  the  phantom.  The  height  of  the 
narrow  rectangle  demarcates  the  volume  averaged  for  B1  quantification,  while  the  dashed  box 
demarcates  the  location  of  T2  evaluation.  The  length  of  the  dashed  box  demarcates  the  slice 
package,  while  the  height  of  the  dashed  box  demarcates  the  diameter  of  the  ROl  for  noise- 
insensitive  T2  quantification,  (b)  T2  variation  across  sequential  imaging  slices  at  RF  field  offset 
between  approximately  -20  and  20%.  The  T2  error  is  the  percentage  standard  deviation  of  T2 
measurements  across  all  slices,  relative  to  the  T2  value  in  the  central  slice  at  minimal  RF  offset. 
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Figure  3:  T2prep  acquisitions  and  comparative  FSE  images  acquired  from  a  Cohort  A  subject 
throughout  the  prostate;  column  (A)  Early  echo  time  T2prep  image;  column  (B)  late  echo  time 
images;  and  column  (C)  FSE  image.  The  T2prep  images  are  windowed  between  0  and  1200  and 
0  and  500  signal  units  respectively.  Similarly,  FSE  images  are  windowed  between  0  and  500 
signal  units.  Intra-prostatic  implanted  gold  fiducial  markers  are  visible  as  focal  hypo-intensities 
in  some  figure  panels  (white  arrows  in  column  (A)). 
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Figure  4:  Representative  T2prep  acquisitions,  acquired  at  apex,  mid  (x2),  and  base  locations  of 
the  prostate  using  an  endo-rectal  coil;  Column  (A)  Early  echo  time  images,  windowed  between 
signal  intensities  of  0  and  1000  (256  grey  levels);  Column  (B)  late  echo  time  images,  windowed 
between  signal  intensities  of  0  and  500  (256  grey  levels);  and  (c)  T2  maps,  windowed  between 
T2  values  of  0  and  200ms  (256  grey  levels).  Three  intra-prostatic  implanted  gold  fiducial 
markers  are  visible  as  focal  hypo-intensities  in  Column  (A),  and  demarcated  with  the  black 
arrows. 
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FSE  T2prep  (TE=88ms)  T2prep  map 


Figure  5:  Representative  images  for  Cohort  A  (top  rows)  and  Cohort  B  (bottom  rows), 
displaying  tumors  as  nodular  hypointensities  in  diagnostic  fast-spin-echo,  T2-weighted  T2prep 
images,  and  T2  maps  created  from  T2prep  images  at  echo  times  of  3  and  88ms.  The  T2  maps 
are  windowed  between  0  and  200ms. 
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